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1 Introduction 

Let Q C H N be a bounded domain with a C 2 -smooth boundary dfl and let u C be 
a subdomain. Write Q = Q x (0, T) with T > and write £ = dfl x (0, T). Consider the 
following parabolic equation: 

dy 

— (x,t) + L Q y(x,t) + e(x,t)y(x,t) = f(x,t), in Q = fix (0,T), ^ ^ 

y(x,t) = 0, on E = ffix (0,T), 

where 

N d •• d 
L Q y(x,t) = - —(a l3 (x)—y(x,t)) + c(x)y(x,t) 

i,j=l 3 1 
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is considered as the system operator. Here and in all that follows, we make the following 
regularity assumptions for the coefficients of L : 

(I) : 

N i 

a«(x) G Lip(U), a^(x) = a^(x), and A*|£| 2 < £ a ij (x)Uj < t;\Z\ 2 , for £ e 

ij=i A 

(1.2) 

with A* a certain positive constant; 

(II) : 

c(x) G L°°(n), e(x,t) G L°°(0, T; L q (Q)) with g > max{iV,2}, and f(x,t) G L 2 (Q). 

(1.3) 

In such a system, we regard e(x, t) as a perturbation in the system conductivity. Suppose 
in the ideal case, namely, in the case when the perturbation e(x,t) = 0, (1.1) has a 
periodic solution y (x,t): 



-t^tOM) + L y (x,t) = f(x,t), in Q, 

y o {x,t)=0, on S, 

y o (x,0)=y o (x,T), in Q. 



;i-4) 



Then the presence of the error term e(x, t) may well destroy the periodicity of the system. 
Indeed, (1.1) may no longer have any periodic solution. (See Section 3.) The problem 
that we are interested in in this paper is to understand if there is a finite (constructible) 
dimensional subspace U C L 2 (Q), such that, after imposing a control u e G U, we can 
restore the periodic solution y e . Moreover, we would like to know if y e is close to y and if 
the energy of u e is small, when e(x, t) is small. Our main purpose of this paper is to show 
that we can indeed achieve this goal in the small perturbation case, even if the control is 
only imposed over a subregion u of Q. The basic tool for this study is the existence and 
energy estimate for the approximate periodic solutions obtained in the author's previous 
paper [8]. 

To state our results, we first recall the definition of approximate periodic solutions 
with respect to the elliptic operator L . 

Notice that L is a symmetric operator. Consider the eigenvalue problem of L : 

f L X(x) = \X(x), 

\X(x)\ Ba = 0. {Lb} 

Making use of the regularity assumptions of the coefficients of L , we know (see, [2] [5], 
for example) that (1.5) has a complete set of eigenvalues {Xj}^ with the associated 
eigenvectors {X j (x)}°^ 1 such that 

L Q Xj{x) = \jXj{x), 



2 



-oo < Ai < A 2 < • • ■ < Xj < ■ ■ ■ < oo, lim Aj = oo, Xj(x) G H^(Q) n C(Q). 

j— >oo 

Choose such that it forms an orthonormal basis of L 2 (Q). Therefore, for any 

oo 

y(x,t) G L 2 (Q), we have y(x,t) = yj(t)Xj(x), where 



3=1 



Vj (t) = (y(x,t),X 3 (x)) = J a y{x,t)Xj{x)dx G L 2 (0,T). 



Definition 1.1. We call y(x,t) is a K- approximate periodic solution of (1.1) with 
respect to L if 

(a): y G C([0, T]; L 2 (fi)) n L 2 (0, T; a weaA; so/nton o/ fi.ij; 

(7>,): y G Sk, where Sk «s £/ie space o/ £/ie following functions: 

S K = {y(x,t)eL 2 (Q); y j (0)=y j (T), forj>K + l, y 3 (t) = J^y(x,t)X,(x)dx}. 



When K = 0, we will always regard X)° =1 = 0. Hence, a 0-approximate periodic 
solution of (1.1) is a regular periodic solution. In what follows, we write (y(-,t),y(-,t)) = 

J y 2 (x,t)dx = t)|| 2 , and we denote y t for the derivative of y(x,t) with respect to t. 
Our first result of this paper can be stated as follows: 

Theorem 1.1 Consider the system (1.1), where e(x,t) is regarded as a perturbation 

in the system conductivity. Suppose that (1.1) has a periodic solution yo(x,t) at the ideal 

case with e(x,t) = 0. Assume that \\e(x,t)\\L°°(o,T;L^(p)) — ess SU P Il e (") OlU^n) < e , 

te(o,T) 

where e < 1 is a small constant which depends only on Lq, Q, N, q, T with q > max{iV, 2}. 
Then there are a non-negative integer K , depending only on L ,Q, N,q,T (but not f), 
and a unique outside force of the form 

K 

u e (x) :^^2u j X j (x) G U = span R {X 1 (x),X 2 (x), ■ ■ ■ ,X Ko (x)}, 
j'=i 

where Uj G R ; such that the following has a unique periodic solution y satisfying: 

( + L y(x,t) + e(x,t)y(x,t) = f(x,t) +u e (x), in Q, 

y(x, t) = 0, on E, 

(y(x,0),^(x)) = (yo^O),^-^)), for j < K , 

k y(x,0) = y(x,T), in Vt. 

(1.6) 
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Moreover, we have the following energy estimate: 



and 



sup \\(y-y )(.,t)\\ 2 + [ \\V(y-y )(;t)\\ 2 dt 
te[o,T] Jo 

< C(system,K )\\e(x,t)\\ 2 LOO{0T . Lq{n)) (l + \a\ 2 + / fdxdt), 
I«e||i2 ( n) < C(system,K )\\e(x,t)\\ 2 L oo {0jT . Lqm (l + \a\ 2 + / fdxdt), (1.8) 

J (J 



where a = (a 1 , a 2 , • • • , a Ko ) = ({y {x, 0), X ± (x)), (y (x, 0), X 2 (x)), ■ ■ ■ , (y (x, 0), X Ko (x))). 
Here and in what follows, C (system, Kq) denotes a constant depending only on Lq, Q, N, q, T , 
which may be different in different contexts. 



In Section 3 of this paper, we will construct an example, showing that without outside 
controls, (1.1) has no periodic solutions in general. This is one of the main features in our 
Theorem 1.1: The control can always be taken from a certain fixed constructible finite 
dimensional subspace to regain the periodicity, while the perturbation space for e(x, t), 
which destroys the periodicity, is of infinite dimension. We also notice that our system 
operator Lq is not assumed to be positive. 

The second part of this work is to consider the same problem as studied in the first part, 
but with the control only imposed over a subregion oo C Q and time interval E C [0, T], 
m(E) > 0. We will similarly obtain the following: 

Theorem 1.2. Suppose that the system (1.1) has a periodic solution yo(x,t) at the 
ideal case with e(x,t) = 0. Then there are a positive integer K , a small constant e > 0, 
depending only on L ,Q, N,q,T (q > max{iV, 2}) 7 such that, when 

\\e(x,t)\\ L oo {0tT . Lg{n)) = ess sup \\e(x,t)\\ Lq(n) < e, 

te(o,T) 

the following has a unique periodic solution: 

dy(x t) K ° 

— — — + L y(x,t) +e(x,t)y(x,t) = f(x,t) +^2xU x )XE(t)u j X j (x), in Q, 

m i=i 
y(x, t) = 0, on S, 

(y(x, Q),Xj(x)) = dj, for j < K , 

y e S Ko , 

(1.9) 

where (a ± , a 2 , • • • , a Ko ) = ((y (x, 0), X^x)), (y (x, 0),X 2 (x)), (y (x, 0),X Ko (x))) = a, 
(ui, u 2 ,---, uk ) = u G R^ . Moreover, 

\u\ 2 <C{system,K,,u)^^^ [ fdxdt), (1.10) 

(m(E)Y Jq 
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and 

SUp \\(y-y )(.,t)\\ 2 + [ T \\V(y-y )(;t)fdt 
* e[ °' T] 2 (1.11) 

< C(system,K ,u) — — K — — — (1 + a + / / dxdt). 

{m{E)Y Jq 

Here, 

Xw(x), XE(t) 

are the characteristic functions for uj and E, respectively; and C '{system, K ,u) is a 
constant depending only on u, L , Q, N, q, T. 

Theorem 1.1 and Theorem 1.2 give stabilization results for the periodic solutions of a 
linear parabolic system under small perturbation of the system conductivity, modifying a 
control from a fixed finite dimensional subspace. We do not know if similar results as in 
Theorem 1.1 hold under the large perturbation case. 

The paper is organized as follows. In Section 2, we prove Theorem 1.1. In Section 3, 
we give an example to show that with a small perturbation e(x,t), (1.1) has no periodic 
solution in general. In section 4, we give the proof of Theorem 1.2. 



2 Small perturbation 

In this Section, we give a proof of Theorem 1.1, based on the author's previous paper 
[8]. For convenience of the reader, we first recall the following result of [8], which will be 
used here. 

Theorem 2.1. Assume (1.2) and (1.3). Let e(x,t) G Ai(q,M), where, for any 
positive number M and q > y, 

M(q,M) := {e(x,t) G L°°(0, T; L q (Q)); ess sup te{0jT) \\e(x,t)\\ Lq{Q) < M}. 

Then, there exists an integer K (L , M, fl, q, N, T) > 0, depending only on (L , M, Q, q, N, T) 
(but not f(x,t)), such that for any K > K (L , M,£l,q, N,T) and any initial value 
a — (ai, a2, • • • , ok) G R k , we have a unique solution to the following equation: 

■ dy(x t) 

— ^ h L y(x,t) + e(x,t)y(x,t) = f(x,t), in Q, 

y(x,t) = 0, on S, (2.1) 

(y(x,0),Xj(x)) = dj, for j < K, 

{ y e S K . 

Moreover, for such a solution y(x,t), we have the following energy estimate: 

sup \\y(,t)\\ 2 + f T \\Vy(,t)\\ 2 dt<C(L ,M,n,q,N,T)(\a\ 2 + f fdxdt). (2.2) 
te[o,r] Jo Jq 
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Now, suppose 2/0 is a periodic solution of (1.1) with e(x,t) = 0, namely, 



-Qf(x, t) + L y (x, t) = f(x, t), 

y (x,t) = 0, 
k Vo(x,Q) = y (x,T), 



in Q, 

on E, 
in Q. 



(2.3) 



Let 



S = Mj(0) = (y (x, 0),Xj(x)), for j = 1, 2, • • • . 

In all that follows, we assume that e(x,t) e .M(g, M) with M = 1. By Theorem 
2.1, there exists an integer K (L , M,il, q, N,T) > 0, such that for the initial value 
a = (ai, a 2 , • • • , a,K ) £ R^ , we have a unique solution y(x,t) satisfying the following 
equations: 



in Q, 

on E, 
for j < K , 



dy<K ^ ~ + £oJ/(z, + e(x, i)y(a:, t) = /(x, t) + ^UjX^x), 

y(x,t) = 0, 

0),Xj(x)) = dj, 

. y e S Ko . 

Here, -u = (ui, tt 2 , • • • , m^ ) G 

Subtracting (2.3) from (2.4), we get the following equation: 

K 

(y ~ Uo)t + L o(y - Do) + e(x,t)(y - y Q ) = J2 u j X j( x ) ~ e(x,t)y , in Q, 
(y(x,t) - yo(x,t)) = 0, 

(y - yo^O) = (y(x, 0) - y (x, 0), Xj(x)) = 0, 
(y-y ) e S Ko . 



(2.4) 



on E, 
for j < K , 



(2.5) 



J : R*° > 



R 



We define a map 
by 

J(ui, u 2 , ■ ■ ■ , u Ko ) = ((y - y )i(T), (y - y ) 2 (T), ■■■ ) {y- y ) Ko (T)). 

Write v = y — y = v + v u . Here, v and v u are the solution of the following equations, 
respectively, 



(v ) t + L v + e(x, t)v = -e(x, t)y , in Q, 

vo = 0, on E, 

(uo)i(0) = 0, for j < K , 
. vo e S Ko - 



(2.6) 
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and 



K 

(v u )t + L v u + e(x, t)v u = ^2ujXj(x), in Q, 



v u = 0, 



K),(o) = o, 

[ v u e S Ko - 



on E, 

for j < K , 



(2.7) 



We are led to the question to find out if there is a vector u = (ui, u 2 , • • • ,%„) G R Ko 
such that 

J(u) = ((y - y )i(n (y - y ) 2 (T), ■ ■ ■ , (y - y ) Ko (T)) = (0, 0, • • • , 0). 

Indeed, if this is the case, then y is a periodic solution with the required estimate as we 
will see later. 

For this purpose, we write J = ((v )i(T), (v ) 2 (T), (v ) Ko (T)) and 



J*(u) = (K)i(T), (v u ) 2 (T), • • • , K)ko(^)). 



Then 



J(u) = J + J*(u). 

Now, it is easy to see that J* is linear in (u\, u 2 ,---, uk )- We next claim that J* is invert- 
ible under the small perturbation case. If not, we can find a vector £ = £ 2) • • • , £k" ) £ 



R x ° with |£| = + £f H 1- £jr = 1 such that J*(£) = 0. Hence, we have a unique 

solution to the following problem: 



K 

w t + L w + e(x, t)w = ^2^jXj(x), in Q, 

j'=i 

w = 0, on S, 

^•(0) = Wj (T) = 0, for j < X , 

s w e S Ko - 



(2.8) 



First, by the energy estimate in Theorem 2.1, we have for w(x,t), 

sup |H-,0H 2 + [ T \\X7w(-,t)\\ 2 dt <C(system)-T-\g\< 
te[o,T] Jo 



< C (system, K ). 



(2.9) 



As mentioned before, we use C (system, K ) to denote a constant depending only on 
Lo, M, Q, q, N, T, which may be different in different contexts. 

00 

Write w = ^2wj(t)Xj(x) as before. Then we have 
3=1 

dw 3(t) 
dt 



+ XjWj(t) + / e(x, t)w(x, t)Xj(x)dx = for j = 1, 2, • • • , K . (2.10) 
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Next, by the Holder inequality (see Claim 2.2 of [8]), we have 



J \e(x,t)w(x,t)X j (x)\dx < C(n,N,q)\\e(x,t)\\ L °o {0iT . Lqm [\\w(-,t)\\ 2 L 2 {n) 

+\\Xj(*)\\ha) + HV^(-,i)||i 2( n) + l|V^.(x)||i 2(n) ] 
< C(Q, N, q) \\e(x, t) |Uoc (0) r;w(n)) [1 + 
+\W;t)\\hw + \\Vw(;t)\\h { n)}- 

By (2.9), we have 

/ / \e(x,t)w(x,t)X j (x)\dxdt < C(system, K )\\e(x,t)\\ L oo^ T . Lg ^))- (2.11) 

v v 

Next, from (2.10), we get 

(e A '*Wj(*))t + / e ( x i t)w(x, t)Xj(x)e Xjt dx = e Aj %-, for j = 1, 2, • • • , 
Integrating the above over [0,T], we get, for j = 1, 2, • • • , K , 

° + F I e(x,t)w(x,t)X j (x)e Xjt dxdt = & Fe Xjt dt. 
Jo Jn Jo 



Namely, 



/ / e(x,t)w(x,t)X j (x)e X]t dxdt 
— ^t_ 1) > for A, ^ 0, 

ne(x, t)w(x, t)Xj(x)dxdt 
i 



T 



for Xj = 0. 



Hence, we get, for j — 1, 2, • • • , if , 



101 < C (system, K ) / |e(x, t)w(x, t)Xj(x)\dxdt 

J •/ 

< C(sj/s*em,ifo)||e(x,*)||Loo( 0> r;Lfl(n))- 



We get 

— # j 

1 = |Cl < C(system,K )^/ir ||e(a;,t)|| L oc (0iT;L9(n)) . 
This gives a contradiction when 

1 



|e(x,i)||Loo( 0j T;L9(n)) < 



C (system, K )^/K^ 



(2.12) 



Therefore, we showed that J* is invertible when ||e(x, i)IU°°(o,T;L«(n)) < e with a certain e 
depending only on L , f2, iV, q, T. 
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Hence, for any given b — (bi, b 2 , ■ ■ ■ , bK ) G R^ , there exists a unique 

u = {u u u 2 , ■ ■ ■ ,u Ko ) e 

such that 

J*(u) = J*(ui,u 2 , ■ ■ ■ ,u Ko ) = (biMr • • ,b Ko )- 
Back to the equation (2.7), we have 

d(v v )j(t) 



Then 



dt 



+ ^j{v u )j(t) + / e(x, t)v u (x, t)Xj(x)dx = Uj, for j = 1, 2, • • • , K . 



dt 



+ J e(x, t)v u (x, t)Xj(x)e X:it dx = Uje Xjt : , for j — 1, 2, • • • , K . 
Integrating the above over [0,T], by the definition of J*, we have 

6 7 e AjT - + / / e(x, t)vJx, t)X j (x)e X]t dxdt = u j [ e Xjt dt, for j = 1, 2, • • • , 
jo Jn JO 

We then get 

bje XjT + [ [ e(x,t)v u (x,t)X j (x)e Xjt dxdt 
Jo Jn 



e(x, t)v u (x, t)Xj(x)dxdt 



for A,- ^ 0, 



for Xj = 0. 



\u 3 \ 2 <2e 2X Ko T \b J \ 2 + 2e 2X Ko T [J^ ^e(x,t)v u (x,t)X j (x)dxdt} 2 

< 2eW|^f + 2e 2X «« T • supjX,| 2 [ / T ||e(, t) || L , (n) K(-, t) \\ Lq , {n) dtf 



Here 1/q + 1/q' = 1. Since f2 is bounded and g' = < 2, by the Holder inequality, we 
have \\v u \\ Lq > {n) < C(tt, q)\\v u \\ L 2 (n) . Hence, 

rT 

[ / ||e(-, t) \\ LH n) \M; t)\\ Lq > (n) dt} 2 < C(n, T, g)||e||^ (0iT;Lq(n)) K|||2 (g) . 

J 

By the energy estimate in Theorem 2.1, we have ||f u ||x,2(Q) < C (system, K )\u\ 2 . 

Hence, as argument before, when ||e|||oo( T\L9(«)) * s sman ) we can solve the above to 
obtain the following: 



\u\ 2 < C(system,K )\b\' 



(2.13) 
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Back to (2.5), we need to find u = (ui, u 2 , - • • , uk ) such that the solution in (2.5) has 
the property (y — y )j(T) = for j = 1, 2, • • • , K . As mentioned before, v = y — yo is 
then a periodic solution. Thus y = v + yo is a periodic solution of (2.4) after applying the 

K 

control force y~]ujXj(x). To this aim, we need only to find u such that 

3=1 

J(u) = or J* (it) = -J . 
By the definition of Jo, Jo = —b = (—b\, — 6 2 , • • • , — &jr ) is given by 

(u )t + J^o + e(x,t)w = -e(x,t)y , in Q, 
f = 0, on S, 

(uo)i(0) = 0, (v )j(T) = -bj, for j < K , 

vo e S Ko - 

By the energy estimate of Theorem 2.1, we have 

\b? <\\vo(;T)\\l Hn) 

< C (system, Kq) \ \ (—ey ) 2 dxdt 

Jo T Jn 

<C(system,K )[ {\\e(x,t)\\ 2 Lq[n)) \\y (-,t)\\ 2 ^ }dt (2.14) 
Jo v " L^(n) 

< C(system,Ko)\\e(x,t)\\l x ^ T . Lqm \\Vyo\\ 2 L 2 (Q) 

< C(system,Ko)\\e(x,t)\\ 2 L o 0( o iT . Lqm (\a\ 2 + / fdxdt), 

where a= (a u a 2 , ■ ■ ■ , a Ko ) = ({y o (x,0),X 1 (x)),{y o (x, 0), X 2 (x)), ■ ■ ■ , (y (x, 0),X Ko (x))). 
Thus, by (2.13), we get 

\u\ 2 < C(system,Ko)\\e(x,t)\\ 2 Loo{0iT . Lqm (l + \a\ 2 + f fdxdt). (2.15) 

J Q 

By (2.2), (2.14) and (2.15), we obtain 

sup \\(y-y )(.,t)\\ 2 + j T \\V(y-yo)(;t)\\ 2 dt 
te[o,T] Jo 

< C(system,Ko)\\e(x,t)\\ 2 LOO{OT . Lqm (l + \a\ 2 + / fdxdt), 
Summarizing the above, we complete the proof of Theorem 1.1. I 
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3 An example 

In this section, we present an example, showing that with a small perturbation e(x,t), 
(1.1) has no periodic solution in general. This demonstrates the importance of an outside 
control to gain back the periodicity as in Theorem 1.1. 

We consider the following one dimensional parabolic equation: 

f Vt ~ Vxx ~ V ~ e{x)y = f(x), < x < vr, < t < T, 

\ y(0,t) = y(7r,t) = 0, 0<t<T. 1 ' > 

Let L e y = —y xx — y — e(x)y with e(x) G C°[0, n]. Suppose is an eigenvalue of L e with 
eigenvectors {Xj(x)} 1 J^ 1 . Then (3.1) has a periodic solution if and only if 



I" f(x)Xj(x)dx = 0, for j = 1,2, 
Jo 



m. 



Now, when e(x) = 0, then is the first eigenvalue of L Q with sinx as a basis of the 
0-eigenspace. Hence, (3.1) has a periodic solution if and only if 

„ n oo oo 

/ f(x) sinxdx = or f(x) = V]aj sinja;, V" \aA 2 < oo. 

J ° 3=2 3=2 

Now suppose e(x) ~ 0. The first eigenvalue A e of L e is given by 



A e = min Je{.Vi V 9 )? 

V eH l J(0,7r),||vp|| L 2 (0 ^ ) = l 



where 



PIT 

J e (<p, if) = / ((fil -Lp 2 - e(x)cp 2 )dx. 
Jo 

(See [5]). Hence, 



A e < min 

>p€Hi(0,n),\\<p\\ L 2 (0tn y 

< + max |e(x)| 

< max |e(x)|. 



/*7T f'7T 

/ (cp 2 — (p 2 )dx + max \e(x)\ / (p 2 dx 
--iJo Jo 



(3.2) 



A e = Je((Pe, <Pe) = / (<Pe)ldx - (1 + e(x))<f%dx 

Jo Jo 

with </? e the eigenvector corresponding to A e and ||</2 e ||L 2 (o,7r) = 1- 
Since is the first eigenvalue of L , we have 



> — max \e(x) 



(3.3) 
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By (3.2) and (3.3), we get 

|A e | < max |e(x)|, and A e — > as e(x) — > 0. 

Next, consider the system with e(x) + A e as the perturbation in the system conductivity: 

f yt-yxx-y-(e(x) + \ e )y = f(x), < x < 7T, < t < T, , . 

\ t/(0,t)=y(7r,t)=0, 0<t<T. 1 • J 

Then when e(x) « 0, we have (e(x) + A e ) « 0. However, if (3.4) still has a periodic 
solution, we have 

PIT 

/ f(x)ip e dx = 0. 

JO 

If this is the case for any given /, we then have 

rn 

/ sin jx(p e dx = 0, for j — 2, 3, • • • . 
Jo 

This implies that ip e = C sin x and thus 

— e(x)smx = A e sina;, or e(x) = — A e . 

This is a contradiction unless e(x) = const.. This shows that for any non-constant small 
perturbation in e(x), for most a priori given /, the periodicity of the system will get lost. 



4 Local stabilization 



In this section, we consider the same problem as studied in Section 2, but with the 
control only imposed over a subregion uj C Q and time interval E C [0, T] with m(E) > 0. 

For the proof of Theorem 1.2, we need the following lemma, whose quantitative version 
in the Laplacian case can be found in [6] and [7]: 



Lemma 4.1 Let Xij(u) = / Xi(x)Xj(x)dx. Then the symmetric matrix X (u, k) = 
(Xjj(w))i<jj<fc is positive definite for any k > 1. In particular, it is invertible. 
Proof of Lemma J^.l: Let a = (ai, a<i, • • ■ , a*,) G R k and let 

r k 

I (a, a) 



Jul j = 1 



dx. 



Then 



I (a, a) = a ■ X(u>, k) • a T , where a T 



( ai \ 

a 2 
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Apparently, I(a,a) > 0. If X(u,k) is not positive definite, then there is a vector a' = 
(a' 1: a' 2) ■ ■ ■ , a' k ) ^ such that I (a', a') = 0. Without loss of generality, assume that a' k ^ 0. 
Hence, 

k 

, £af j X j (x)\ u = 0. (4.1) 
i=i 

We thus get over u: 



a' 



X k (x) = ]T bjXjix), with b 3 = - -f . (4.2) 

i<fc 

Applying (L ) m to (4.2) over a;, we have 



j<k a k 



Kx k { x ) = Y,h>yxAx)- 

j<k 



We get 

X k (x) = 

j<k 



X k (x) = '£b j &) m X j (x) over u. 



X k 



Letting m — > oo, we get over u 



X k (x)= £ bjXjix), (4.3) 

k'<j<k 



where 

f Aj = A fc , for j > fc', 
| Xj < A fc , for j < k'. 

By (4.4), we get over Q, 



(4.4) 



L (Xfc(x)- J2 b 3 X 3 (x)) = \ k X k (x) - J2 bjXjXjix) = X k [X k (x) - £ b 3 X 3 (x)]. 

k'<j<k k'<j<l k'<j<k 

By (4.3) and the unique continuation for solutions of elliptic equations, we get 

X k (x) — ^ bjX 3 (x) = over f2. 

k'<j<k 

This contradicts the linear independence of the system {Xj}M 

Proof of Theorem 1.2.: Similar to the proof of Theorem 1.2, we need only to find 
a vector u = (u±, u 2 , ■ ■ ■ , uk ) G R^ such that 

where 

^C(^) — ((^j M 3 ^ ^(z)), • • • , {v(x, T),X Ko (x))) — ('Ui(T), v 2 (T), • • • , v Ko (T)) 
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with v the solution of the following equation: 

v t + L v + e(x,t)v = Y,xU x )XE(t)u j X j (x), in Q, 



3=1 



v = 0, 
Vj(0) = 0, 
v e S Ko - 



on £, 

for j <K , 



(4.5) 



and 



\w = ((v )i(T), (v ) 2 (T), • • • , (v )k (T)) 
with vq the solution of the following system 

(v ) t + L v + e(x, t)v = -e(x, t)y , in Q, 

Vq = 0, on E, 

(^■(0) = 0, for j < K , 
. vo e S Ko - 



(4.6) 



— * — * 

In the same way, if J* is not invertible, then for a vector £ = £ 2 , • • • , £k ) with |£| = 1, 
we have a solution to the following system: 

K 

v t + L v + e(x, t)v = ^xMXE^jX^x), in Q, 

3=1 

v = 0, on S, 

v j (0) = = v j (T), forj<K , 
k w e S Ko . 



We then get 

+ V>i(*) + / e(x,t)v(x,t)X,-(x)dx = 5Z^Xb(*)-^«M> for j = l,2,---,iT - 
We similarly get 

(e x ^(t))' t + e A ^ ^ e(x, t)v(x, t)Xj{x)dx = ]T 6x B (f)X,», for j = 1, 2, • • • , tf . 



+ 

We then get 

I f T ^ 



e Xjt e(x,t)v(x,t)X j (x)dxdt = J e Xjt J2jiXE(t)Xij(u;)dt. 



I e x ^ X E(t)dt 
Jo 



f e X2t XE (t)dt 
Jo 



j T e x ^x E {t)dt 
Jo 



X(u,K ) 



( 6 \ 

v u, y 
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o T Jn 



e 1 e(x, t)v(x, i)Xi(x)dxdt 



f f e X2t e(x,t)v(x,t)X 2 (x)dxdt 
Jo Jn 



o Jn 



e XK ^e(x, t)v(x, t)X Ko (x)dxdt 



I 



( 6 \ 

v 6r y 



= x(^,k )- 1 



e Xlt XE (t)dt)- 1 f e x ^evX x dxdt ' 

e X2t x E {t)dty l ( e X2t evX 2 dxdt 
JO. 



\ 



(4.7) 



/ 



e XK » l x E {t)dty l f e Xli ^evX Ko dxdt 
io Jq 

By Lemma 4.1, we know X(u>, Kq)~ 1 is a bounded linear operator from H K " to R^ . 
By the energy estimate in Theorem 2.1, we have for v(x,t), 

sup \\v{-,t)\\ 2 + / \\Vv{-,t)\\ 2 dt < C{system, K ) / Q2x u (x)x E (t)^X j (x)) 2 dxdt 
te[o,T] Jo Jq j=1 

<C(system,K )T\Z\ 2 
< C (system, K ). 

By the Holder inequality, we have, 

J \evXj\dx< C(fi,iV,g)||e|| L o O(0iT . L9(f , )) [l + A J 2 + \\v(-,t)\\ 2 + ||Vu(-,f)|| : 
Together with (4.8), we thus have 

/ / \evXj\dxdt < C(system,K )\\e\\ L oo^ iT . Lq{n )y 

« « f2 

Back to (4.7), we have 

|£| 2 < ^(system,^,^)^^!!^^,^)^ 1 !! 2 ^!!!^^.^^)) 
< C(system, u, K ) ^ E ^ 2 \\ e \\L™{o,T;Li(n))- 

Hence, when Hellioo^T^rn)) ^ s sufficient small, we get \£\ 2 < 1. This gives a contradic- 
tion. Therefore, we showed that J* is invertible under small perturbation. By the same 
arguments as those in the proof of Theorem 1.1, we can also show the energy estimates 
as stated in Theorem 1.2. This completes the proof of Theorem 1.2. ■ 



(4.8) 
(4.9) 

(4.10) 



15 



References 



[1] V. Barbu and G. S. Wang, Feedback stabilization of periodic solutions to nonlinear 
parabolic-like evolution systems. Indiana Univ. Math. J. 54 (2005), no. 6, 1521-1546. 

[2] L. C. Evans, Partial Differential Equations, Graduate Studies in Mathematics (Vol 
19), American Mathematical Society, 1998. 

[3] A. Friedman, S. Y. Huang and J. M. Yong, Optimal periodic control for the two-phase 
Stefan problem. SIAM J. Control Optim. 26 (1988), no. 1, 23-41. 

[4] A. Friedman and L. S. Jiang, Periodic solutions for a thermostat control problem. 
Comm. Partial Differential Equations 13 (1988), no. 5, 515-550. 

[5] D. Gilbarg and N. S.Trudinger, Elliptic Partial Differential Equations of Second 
Order. Grundlehren der mathematischen Wissenschaften, a Series of Comprehensive 
Studies in Mathematics, Vol. 224, 2. Edition, Springer 1977, 1983. 

[6] G. Lebeau and L. Robbiano, Contrdle exact de l'equation de la chaleur. Comm. PDE, 
20(1-2), 335-356, 1995. 

[7] G. Lebeau and E. Zuazua, Null-controllability of a system of linear thermoelasticity. 
Arch. Rational Mech.Anal. 141, 297-329, 1998. 

[8] L. Lei, Identification of Parameters through the Approximate Periodic Solutions of a 
Parabolic System, Journal of Optimization Theory and Applications, Vol. 137, PP: 
185-204, 2008. 



16 



